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Abstract 

A cosmological pseudoscalar field coupled to hypercharge topological num- 
ber density can exponentially amplify hyperelectric and hypermagnetic fields 
while coherently rolling or oscillating, leading to the formation of a time- 
dependent condensate of topological number density. The topological con- 
densate can be converted, under certain conditions, into baryons in sufficient 
quantity to explain the observed baryon asymmetry in the universe. The am- 
plified hypermagnetic field can perhaps sufficiently strengthen the electroweak 
phase transition, and by doing so, save any pre-existing baryon number asym- 
metry from extinction. 
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To generate baryon asymmetry from a state in which quark and antiquarks have equal 
abundances, the three celebrated Sakharov conditions have to be satisfied. In the scenario 
we propose, the Sakharov conditions are satisfied during an epoch in which a cosmological 
pseudo-scalar field coupled to hypercharge topological number density coherently rolls or 
oscillates. The scalar field motion creates a time dependent hypercharge topological num- 
ber condensate which violates fermion number conservation through the abelian anomalous 
coupling PJ, and establishes the required departure from equilibrium, and its interactions 
violate C and CP symmetries. Pseudoscalar fields with the proposed axion-like coupling 
appear in several possible extensions of the Standard Model. They typically have only 
perturbative derivative interactions and therefore vanishing potential at high temperatures, 
and acquire a potential at lower temperatures through non-perturbative interactions. Their 
potentials take the generic form V(4>) = VqV(4>/ /), where V is a bounded periodic function 
characterized by a scale / ("Peccei-Quinn" scale), which could be as high as the Planck 
scale and a much smaller mass m ~ V 2 / /, which could be as low as a fraction of an eV, or 
as high as 10 12 GeV. A particularly interesting mass range is the TeV range, expected to 
appear if potential generation is associated with supersymmetry breaking. 

The fundamental role of hypermagnetic and hyperelectric fields in the context of elec- 
troweak (EW) baryogenesis has been recognized recently in several investigations Rflj.It was 
observed that: i) a topological number condensate can be released at the EW transition in 
the form of leptons and baryons, and ii) strong enough hypermagnetic fields could make the 
EW transition strongly first order even for those large values of the Higgs mass that have 
not been ruled out by LEP II experiment. Scalars with axion-like coupling to hypercharge 
fields were previously considered in Amplification of ordinary electromagnetic (EM) 

fields by such scalar fields was discussed in j| and their possible use for baryogenesis in ||. 

We will assume that the universe is homogeneous and isotropic, and can be described by 
a conformally flat metric ds 2 = a 2 (r])(dr] 2 — dx\ — dx\ — dxf) where a(rj) is the scale factor of 
the universe, and rj is conformal time related to cosmic time t as a(7])di] = dt. In addition 
to the standard model fields we will consider a time-dependent pseudoscalar field 4>(rj) with 
coupling \4>YY to the U(l)y hypercharge field strength and a potential V{4>) = VqV(<P/ f). 
The coupling constant A (which we will take as positive) has units of mass -1 A ~ 1/M. 
For QCD axions, M ~ /, but in general, it is not always the case, and we will therefore 
take M to be a free parameter, and in particular allow M < f. We will assume that the 
universe is radiation dominated at some early time rj = at T > lOOGeV^ before the scalar 
field becomes relevant. 

Maxwell's equations describing the hyper EM fields (we will drop the hyper from now 
on for brevity) in the resistive approximation |3j of the EW plasma, coupled to the heavy 
pseudoscalar are the following 

(z) V • E = (ii) V • B = 

-» -> dB -> 

(izi) J = aE (iv) — — = —V x E 

or] 

We have rescaled the electric and magnetic fields E = a 2 (rj)S, B = a 2 (rj)B, and the physical 
conductivity o = a(rj)a c . In the EW plasma a c ~ 10T 0. The fields £, B are the flat 
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space EM fields. We have assumed for simplicity vanishing bulk velocity v of the plasma 
and vanishing chemical potentials for all species. Our results can be easily generalized to 
include non-zero bulk velocity or right electron chemical potential. The equation for the 
pseudoscalar is the following 

d 2 d> dd> 9 dV(d>) , 9 - -> 

y + 2aE-j- + a 2 = Xa 2 E ■ B, (2) 



drj 2 drj 



where H = 4^ is the Hubble parameter. We will neglect the backreaction of the electro- 
magnetic fields on the scalar field since it is irrelevant for most of the physics we would like 
to explore. The cosmic friction term 2aH^ will not be relevant for most of our discussion. 
We therefore solve eq.@ with vanishing r.h.s., and substitute the resulting <p(rj) into eq.(|]). 

We are interested in solutions of the form E(x,t)) = J d 3 k e~ tkx eg 6^(77), B{x,rj) = 
J d 3 k e~ tkx bz /%(?7), f° r which the electric and magnetic modes are parallel to each other. 
We find that the Fourier modes eg and b^ are related, 6g ± = b k (ei ± ie^), ke^e^irj) = 

±b^ ± ^- , where <§i, e 2 are unit vectors in the plane perpendicular to k such that (ex, e 2 , k) is 

a right-handed system. This type of solution has oc Vx^oc and V x E% oc E%. 
The function f3r(r)) obeys the following equation, 

d 2 ^ dpk* ( 2 dep \ ± 

Before plunging into a detailed discussion of the solutions of eq.flU), it is useful to consider 
the simple case of a constant f. Then the solutions are simply f3 ± {j]) = flf e^i ^+13^ e U2 ^ , 
where u;^ 2 are the two roots of the quadratic equation u 2 + au; + {k 2 ± Xfjkj = 0, uf^ = 



-a ± ^a 2 - 4 (p ± X^k 

Now the qualitative behaviour of solutions is clear. If one of the eigenvalues is positive, 
which can happen only if A If- > k, EM fields can grow exponentially. Otherwise fields are 
either oscillating or damped, as in ordinary magnetohydrodynamics. To obtain significant 
amplification, coherent scalar field velocities over a duration are necessary, larger velocities 
leading to larger amplification. We estimate the wavenumber k of maximal amplification by 
looking for the maximum of 10 as a function of k, k max = |A f . 

Another interesting approximate solution relevant for rolling fields can be obtained 
for k ^ T. In this case, eq.(^) can be approximated by a first order equation cr|^ + 

(k 2 ± X^k) p±{rj) = 0, which can be solved exactly, (3 = (3 e - ^ + AA< ^ k / a , where 
A0(t?) = |0(r/) — 0(0) I . The amplified mode is determined by the sign of 0(r/) — 0(0). The 
amplification factor A ± (k, 77) = f3 ± {rf) / (3 ± {Qi) is maximal for k max r] = |AA0, A (k max , 77) = 
e x(^<t>) ^ _ Looking at 77 ~ t]EW we obtain v J wCT ~ 10~ 16 , and therefore to obtain amplifi- 
cation AA0 > 10 8 . For specific models an upper bound on AA0 may appear, narrowing the 
range of allowed parameter space. A value of AA0 > 10 8 is not unnatural, for example, such 
a value is obtained if the typical scale for scalar field motion is the Planck scale, as happens 
in many models of supergravity, and A < l/10 w GeV . 
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FIGURES 
Time Dependence: Magnetic Field Amplification 
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FIG. 1. Amplification of EM fields. The function (3, exponentially scaled, is shown for the 
parameters A = 55, k/m = 20, m = 6TeV and a = 40T at T = ITeV. 

If the scalar field is oscillating, the field's velocity changes sign periodically and both 
modes can be amplified, each during a different part of the cycle. Each mode is amplified 
during one part of the cycle and damped during the other part of the cycle. Net amplification 
results when amplification overcomes damping, which occurs roughly when A^fc > a 2 (recall 
that a ~ 10T). Total amplification is exponential in the number of cycles. 

We have studied numerically the solutions of eq. (|]) for the case of an oscillating field for 
different parameters and situations. The concrete potential that we used for the numerical 
analysis is V(<fi) = |m 2 2 , assuming that the initial amplitude of the field is of order / 
with A = //M > 1, f/m>l. We find results in agreement with the previous qualitative 
discussion. Amplification occurs for a limited range of Fourier modes, peaked around k/m ~ 
A. The modes of the EM fields are oscillating with (sometimes complicated) periodic time 
dependence and an exponentially growing amplitude. In Fig. [l] we show an example of the 
time dependence of amplified EM fields for a specific mode. For the range of parameters in 
which fields are amplified, the amount of amplification per cycle for each of the two modes 
is very well approximated by the same constant T(k/m, A, a). A good approximate estimate 
for the average amplification after N cycles is therefore A ± (k, r/) = Af ± k^ Nr , where Af ± k 
represents the transient influence of the initial conditions of EM and scalar fields. 

The scalar field is a very efficient amplifier of EM fields. For example, to obtain an 
amplification of 10 12 for A = 55, k/m = 20, m = 6TeV and a = 40T, for oscillations 
occurring at a temperature of ITeV, we need just one cycle! In Fig. |2| we have presented an 
example of the dependence of the amplification per cycle factor T, on different parameters. 
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Amplification Rate per Oscillation: 

Dependence on Temperature 
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Amplification Rate per Oscillation: 
Dependence on Scalar Field Amplitude 
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FIG. 2. Amplification per cycle of EM fields, as a function of temperature, wave number, and 
scalar field amplitude. 

A detailed discussion of the end of oscillations or rolling is beyond the scope of the 
current investigation. However, we do know that once the oscillations or rolling stop, the 
fields are no longer amplified and obey a diffusion equation 0. Modes with wave number 
below the diffusion value k < k a T ~ 10~ 8 T, where !— = 1, remain almost constant until 
the EW transition, their amplitude goes down as T 2 , and energy density as T 4 , maintaining 
a constant ratio with the environment radiation. Modes with k/T > k a decay quickly, 
washing out the results of amplification. We have seen that the range of amplified momenta 
for oscillating fields is not too different than T, therefore scalar field oscillations have to 
occur just before, or during the EW transition. In that case, the amplified fields do not 
have enough time to be damped by diffusion. If the field is rolling, momenta fc<T can be 
amplified, and therefore the rolling can end sometime before the EW transition. 

To obtain the average magnetic energy density in amplified fields we have to average 
over the magnetic and scalar fields initial conditions. The averaging has different reasons: i) 
statistical fluctuations inside the horizon of magnetic field initial conditions due to thermal 
and quantum fluctuations, and ii) variations over causally disconnected regions. We will 
assume translation and rotation invariance of initial conditions, (B ± (k,r])B ± (£,r]) }\ v =o — 
5 3 (k — £) 5 + ^f ± (k). Using the magnetic power spectrum Pg(k) = Airk 3 f ± (k), the average 
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initial energy density ps(0), given by ps(0) = < B 2 >|^=o, can be expressed as Pb(0) = 

oo . N 

±fd\nk(P+(k) + P B (k)). 

The two modes of the magnetic field are amplified by amplification factors A^(k, rj), 
which depend on the parameters of the model as discussed previously. The amplification 
factors depend, in addition, on the initial conditions of the magnetic field and on the initial 
conditions of the scalar field. For example, if the velocity of the scalar field d<f)/dn is constant, 
the sign of the velocity determines that only A + or A~ will be non-vanishing. 

The magnetic energy density in amplified fields is given by 



P B = ]tJ dlnk { P B~( k ) \A-(k, v )\ 2 + P+(k)\A + (k, v ) 



(4) 



The Chern-Simons number density is given by Ancs 



2 12 V 

/ drj{E ■ B)., where y R 
o 



16vr 2 



is the hypercharge of the right electron, and g' is the hypercharge gauge coupling. Since 
in our case Ef(rj) = ±\d v B${rj), then (E ± (k,r ] )B ± (£,r ] )*) = ±Hd r ,(B ± (k,rj)B ± (£,rj)*). 
Neglecting the initial ncs density, 



ncs 



via' 2 

327T 2 



J d\nk^P B (k) \A-(k, V )\ 2 -P^(k)\A + (k,v)\ 2 Y 



(5) 



Let us define two k-dependent asymmetry parameters, 7^ 5 



and jl s 



p s( fc )+ p K fc )' 

|^-(M)^+)^4+(fc'^)| a ' Both parameters vary from - no asymmetry, to ±1 - maximum asym- 



metry. Using the asymmetry parameters we may relate ncs to 

n cs ^Jffdlnkl PB{ k)^f±^. 
4vr J k I + Castas 



(6) 



Further, we can compute the fractional energy density in coherent magnetic field configu- 
rations VlB(k) = pB(k)/p c , where p c is the critical energy density, and the Chern-Simons 
fractional number density ncs/s, where s is the entropy density, and compare them. If the 

§q<7*T 4 and s 



universe is radiation dominated p, 

yh' 2 H 



45 



■g* s T 3 , therefore 



ncs 



4vr Ag, 



d^B(k)^pf. 
k 1 + IasTas 



(7) 



If the amplification factors A k are, as we have seen, sharp functions of k, peaked at k Ti 
and have width in k of approximately k max , then 

T 



ncs 



B v 
lAS + lAS 

1 + llslts 



0.01 



kr, 



-Q B (k n 



( kmax ) j 



(8) 



where we have used numerical values for the prefactor coefficients. 

According to Giovannini and Shaposhnikov ||, this Chern-Simons number will be re- 
leased in the form of fermions which will not be erased if the EW transition is strongly first 
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order ||, and will generate a baryon asymmetry, — = — f^j 2 -. An equal lepton number 
would also be generated by the same mechanism so that B — L is conserved. Note that the 
fact that baryon number asymmetry is generated at k max ^ does not mean that baryon 
density is actually inhomogeneous on this short length scale L max ~ l/k max . Comoving 
neutron diffusion distance at the beginning of nucleosynthesis is much longer than L max 
PJlO|, so that by that time inhomogeneities would have been erased by free streaming |J. 

If T j k max is not too different than unity, as we have seen for the case of oscillating field, 
and 7^4 S and 7^ 5 are small, it is possible to obtain ^ ~ 10~ 10 and have strong magnetic 
fields Qb ~ 1 present during the EW transition. If T/k max is large and 7^ 5 is order unity 
as we have seen in the rolling case, it is not possible to have strong magnetic fields without 
producing too many baryons. 

The existence of hypermagnetic fields during the EW transition can influence the nature 
of the EW transition. For example, a homogeneous hypermagnetic field adds a pressure 
term to the symmetric phase which could lower the transition temperature and strengthen 
the EW transition. This well known effect in conductor-superconductor phase transitions 
| l"T| was discussed in several investigations @,[|, with seemingly inconclusive results, for the 
moment. 

Our conclusions are therefore that 

1. Depending on the ratio T/k max and on the asymmetry parameters, it is possible to create 
the desired ratio ^j 2 - ~ 10~ 10 , and therefore to generate the observed baryon asymmetry in 
the universe. 

2. Some asymmetry in the initial conditions of either B or is required. The asymmetry 
can be very small and can be induced at a much higher scale. Possible sources for the 
initial asymmetry are temperature dependent potential that traps at a preferred position, 
asymmetry in quantum fluctuations, etc. Large asymmetry in <fi can be expected and appears 
naturally. 
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